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Abstract 

The Kirchhoff 's theory for thin, inextensible, elastic rods with nonhomogeneous cross 
section is studied. The Young's and shear moduli of the rod are considered to vary radially, 
and it is shown that an analytical solution for the constitutive relations can be obtained for 
circular cross section and constant Poisson's ratio. We comment on possible applications 
of our results. 
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Filamentous structures are of great interest in both academic and industrial fields. 
From biological fibers and biopolymers to different kinds of wires and cables, the 
study of statics and dynamics of these systems has brought significant contribu- 
tion to both the understanding of the life and nature, and the development of new 
technological devices. 



The Kirchhoff rod model has been consi dered a good framework to s 



ics and dynamics of filaments in Biology ( Schlick , 



1995 



Olson, 



t udy both stat- 



1996 



Wolgemuth et al. 
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2000; 



Gorielv and Tabor 



Fons eca and de Aguiai 



1998 



Mc Millen and Go rielv 



2002; 



Tobias et al. 



2001) and in Engineering ( Su n and Leonan 



3 mil 



2000; 



Gottlieb ^dPerkins, 



1999). In most cases, the rod is considered to be completely homogeneous. But 



some nonhomogeneities along the rod have been analyzed in the literature. Tridi- 
mensional conformations of nonhomogeneous rods m ay present chaotic behavior 



(Mielke and Holmes, 



1988 



Davies and Moon, 



1993) and deviations from the he- 



lix pattern were sh own to occur in the case of a rod with periodic variation of its 



Young's modulus (Ida Fo nseca et al. 



stiffness nonhomogeneity, 



200 3J). In order to in vestigate the effects of 



da Fonseca and de Aguiarl ([2003) developed a method 



for finding equilibrium solutions of the static Kirchhoff equations for rods sub- 
jected to given boundary conditi ons. A geometric study of ro ds with varying cross 



section radius was performed by 



da Fonseca and Malta! d2003l) . In these examples, 



the rod nonhomogeneities are function only of the arc-length along the rod axis 
(within a section transversal to the axis the system is homogeneous). 



Here we shall deal with nonhomogeneities that are function only of the distance to 
the rod axis. In this case, points lying at a given distance from the rod axis will have 
the same mechanical properties even if they belong to different cross sections. The 
elasticity problem is defined as the specification of the so-calle d state of str ess. It 



requrres 



the knowledge of the stre ss at e very point of the b ody (ILove , 



work by Zha ng and Hasebd (119991) and by 



Chen et al. 



1934). The 



( 2000J) ar e examples of appli- 



cation of Elasticity theory to nonhomogeneous cylindrical rods. lZhang and Hasebe 
( 1999) have obtained an exact solution for the stress of a radially nonhomogeneous 
hollow circular cylinder with exponential radial variation of the Young's modulus 
and constant Poisson's ratio. Here, we address the problem of obtaining an analyt- 
ical solution for the stress of a nonhomogeneous rod within the approximations of 
the Kirchhoff s theory for thin rods. We shall show that an analytical solution can 



2 



be obtained for a rod of circular cross section presenting any kind of radial vari- 
ation of the Young's modulus, but constant Poisson's ratio. Since the Kirchhoff's 
rod theory is largely used in modeling long, thin and inextensible elastic rods, our 
solution can be of interest in a large range of applications. 



Example s of real system s presenting nonhomog eneous cross se ction are coaxial 



cables (ITang et al. 



2001 



stranded DNA molecule (Calladine and Drew 



), coated optical fibers ihietal. . 

1 



2002) and the double 



1999). The process of coating struc- 



tures, where thin layers of a given material are deposited on a given body, is largely 
used in industries as, for example, the process called nitriding t hat improves the so - 



1998). 



called tribo-mechanical properties of engineering components (Miola et al. 
Coating p rocesses are also used in basic research. For example, ISalvadori et al. 



(2003a b) measured the Young's modulus of gold thin films deposited in Atomic 
Force Microscopy (AFM) cantilevers. They found that the Young's modulus of the 
gold thin films is about 12% smaller than its bulk elastic modulus. This information 
is interesting for the analysis we shall present below. 



The results of the Kirchhoff model were recently derived by Mora and Miiller ( 2003) 
through the rigorous method of T-convergence. They showed that the nonlinear 
bending-torsion theory for inextensible rods arises as the T-limit of three-dimensi- 



onal nonlinear theory of elasticity. Nevertheless, we s 



Dill 



lall follow the more simple 



( 1992) since our final solution can 



derivation of the Kirchhoff equations due to 
be directly obtained using this approach. Our aim is to obtain an analytical solution 
for the so-called constitutive relations that relate the components of the moment to 
the components of a vector representing the deformations of the rod. 



Dill (1992) presented the derivation of the Kirchhoff equations from the classical 



conservation laws of linear and angular momentum for a three-dimensional body 
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with a surface area A enclosing a volume V: 

J A P n dS + JyidV = J vP XdV, 

(1) 

J A (X x Pn )dS + J v (Xx f)dV = J v p(X x X)dV , 

where p n is the contact force per unit area exerted on the oriented surface element 
dS = ndS, p is the mass density and X is the position with respect to a fixed 
origin (a dot indicates time derivative). External forces per unit volume acting on 
the body are represented by f . We shall drop the derivatives with respect to time 
since dynamical problems will not be considered here. 

The axis of the rod is defined as a smooth curve x, in the 3D space, parametrized 
by the arc-length s: x = x(s). A director basis is defined at each point of the curve, 
with d 3 chosen as the tangent vector, d 3 = x' (the prime denotes differentiation 
with respect to s). The orthonormal vectors, di and d 2 , lie in the plane normal to 
d 3 . We choose these vectors in such a way that d!,d 2 ,d 3 form a right-handed or- 
thonormal basis at each point of the rod axis. The space variation of the director 
basis along the curve x is controlled by the twist equation: d- = k x dj. The com- 
ponents in the director basis, ki and k 2 , of the twist vector, k, are the components 
of the curvature of the rod, and k 3 is the rod twist density. 

In the Kirchhoff's theory the rod is seen as an assembly of short segments. Each 
segment is loaded by contact forces from the adjacent segments. The equations (1) 
are applied to each segment in order to obtain a one dimensional set of differen- 
tial equations for the static and dynamics of the rod. Each segment of infinitesimal 
thickness will be referred to as a "cross section". The theory assumes that the cross 
section radius, at all points, is much smaller than both the total length and the cur- 
vature of the rod. Since equation (1) involves integration over the volume V of the 
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body or over the area A enclosing the volume V of the body, the nonhomogeneities 
in the cross section do not influence the main Kirchhoff equations. Only the con- 
stitutive relations may depend on the nonhomogeneities in the cross sections. The 



genera 



(ILove . 



equation 



1934 



r elating stress and strain for a homogeneous isotropic material is 



Dill, 



1992) 



S = 2^E + A(Tr£)l 



(2) 



where the tilde is used to denote a tensor. S is the stress tensor, E is the strain 
tensor and 1 is the unitary matrix, u is the shear modulus and A is one of the elastic 



constants of Lame dLove , 



1934; 



Dill, 



19921) . The equation (2) is valid for relatively 



small stresses where the linear theory of elasticity holds true. 



Our problem consists of obtaining the constitutive relations for a rod made of a 
nonhomogeneous isotropic material. In order to clarify the idea of nonhomogeneity 
and isotropy of the material, we must remember that even though the length and 
the radius of curvature of the rod have to be much bigger than the cross section 
radius (for the validity of the Kirchhoff model), the rod is a three-dimensional body 
formed by elements of volume that are considered to be a continuous isotropic 
medium. Nevertheless, some of the elastic properties of the rod can vary from one 
element of volume to another, as for example, in the radial direction. Elasticity 
theory fo r nonhomogeneou s isotropic media has been s tudied since 1960s (see, for 



example (Ro stovtsev 



1964 



Plevako 



1971 



Chen et al. 



2001)) 



Due to the fact that only rods with circular cross sections lead to an analytical 



solution for the constitutive relations, so we shall consider on 



for the elastic properties of the rod. As in dRostovtsev 



1964; 



V radial variation s 



Chen et al. . 



2001), 



the equation (2) remains valid for nonhomogeneous [i and A. Here, fi = /x(r) and 
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A = A(r), where r is the distance of the point to the axis of the rod. 

The point where the axis intersects the plane of the cross section is the origin of a 
Cartesian basis lying in this plane. This Cartesian basis can be defined by the two 
vectors d x and d 2 of the director basis. The Young's modulus E is connected to the 
shear modulus /j, through the Pois son's ratio, v\ v + 1 = E/2/j,. 



The components (i, j = 1,2,3) of th e strain ten sor, E (which is symmetric), 



within the Kirchhoff 's theory, is given by (Dill, 



1992) 



2V,9X (3 " r dX a > ' 



(3) 



e 33 = (*i - kf ] )X 2 - {k 2 - k^X, 



where greek labels equal 1, 2, and Xi and X 2 are the components of the position 
vector of a material point of the rod in the Cartesian basis (di , d 2 ). A;- ^ is the 
i-th component of the twist vector, in the director basis, for the undeformed rod, 
also known as intrinsic curvature. Ui, i = 1,2, 3, are the components of the dis- 
placement of the material point. The solution for = Ui(Xi, X 2 ) constitutes the 
solution for the elasticity problem. 



In order to find the solutions for Ui (i = 1,2, 3) w e shall use the local balance of 



momentum for the components of the stress tensor (Dill, 



1992) 



E 



dS, 



al 



{0X o 



0, Z = 1,2,3 , 



(4) 



where is the ij-th component of the stress tensor, S. The boundary conditions 
are provided by the load conditions on the rod lateral surface. In the Kirchhoff's 
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theory, it is given by 

2 

Y,n a S al = 0, 1 = 1,2,3, (5) 

a=l 

where n x and n 2 are the components of the unit outward vector, normal to the 
boundary of the undeformed cross section. Therefore, the equation (4) must be 
solved for ui subjected to the boundary conditions defined in (5). 

The equation (4) can be separated and solved in two sets. The first one includes all 
terms with the index I = 3 (equation (4)): 

EdS a3 

= ' (6) 

Q=l aX CX 

S a3 = 2/i(r)e a3 , (7) 

e a s = l(§^ + (h-kf)X,). (8) 

The boundary conditions for this set are: 



£ n a S a3 = . (9) 



The distance to the origin, r, is connected to Xi and X 2 through the polar transfor- 



mation X 1 = rcosO, X 2 = r sin 9 (r = yXf + X%). Substituting equations (7) 
and (8) in equation (6), and using the following relation 



Or X a 



dX a 



(10) 



equation (6) reads 



9 ( ( \ du *\ , 9 ( ( \ 9u *\ n ni\ 



dX^ 'dX x ' dX 2 ^ y J dX 2 



1 



This equation describes the torsion of a rod subjected to the 



given by equation (9). The known solution is (|Love . 



u 3 = (k 3 -k$ s >)< P {x 1 ,x 2 ) 



boundary condition 



Dill, 



1992) 



(12) 



where <p(Xi, X 2 ) is known as warping (Dili 1 19921) or torsion fLoveU 19341) func- 
tion, and must satisfy: 



d 2 (p d 2 <p ldfi dip dip 
^dXl + dXl ) + rTr^dXl + X2 dX 2 ] = ° ' 



(13) 



The boundary condition (9) becomes: 

n 1 (^-X 2 )+n 2 (^ + X 1 ) = 



dX 1 



"dX 2 



(14) 



and must be satisfied for all X x and X 2 such that JX 2 + Xf = h, h being the cross 
section radius. 



In the homogeneous case, the function ip depends only on the geometry of the cross 
section. In the nonhomogenoeus case, it also depends on how the shear modulus 
\i varies with r. Nevertheless, by inspection of the equation (14) we see that the 
solution for the homogeneous case with circular cross section (cp(Xi, X 2 ) = 0) also 
satisfies the equation (13) and the boundary condition (14). Therefore, we consider 
here only rods with circular cross section. 



The second set consists of the remaining equations: 
. 



1-y 0X a 



S a p = 2 / u(r)e Q/3 + A(r)(^ e mm )5 al3 , , 



m=l 



1 ( duq dup 
tap- 2 ^dXp dXj ' 



(15) 

(16) 
(17) 
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where /3 — 1, 2 and 5 Q/3 is the Kronecker delta. The boundary conditions for this 
set of equations are: 



(18) 

Eo=l n aS a 2 = , 

If the Poisson's ratio, v, is a constant within the cross section, it is possible to 
show that the solution for u\ and u 2 in eqs. (15-17) have the same form of the 
homogeneous case, even for \x = fi(r) and A = A(r), and the boundary conditions 
(5) are satisfied. Since the stresses considered here are such that the rod remains 
inextensible, the assumption of constant Poisson's ratio is reasonable in spite of 
having varying Young's and shear moduli. In this case, the ratio of Young's over 
shear moduli must satisfy ^j^y = Constant = v + 1. 

The explicit solutions for ui and u 2 are 

u 1 (X 1 ,X 2 ) = -v(ki - M 0) )^2 + - k& 0) )(X* - XI) , 

(19) 

u 2 (X 1 ,X 2 ) = v{k 2 - kf)X x X 2 + \{k x - kf){X\ - XI) . 



Now, we can calculate the constitutive relations in terms of the components of the 
twist vector. The definition of the total moment of the cross section, M, is 



JrxpsdS, (20) 



M = 

s 



where S is the area of the cross section, r is the position vector in the plane of the 
cross section, given by 

r = X 1 d 1 +X 2 d 2 , (21) 
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and P5 is the contact force per unit area in the cross section that, in terms of the 
stress tensor, is given by: 

p s = d 3 .S = S 31 di + S 32 d 2 + S 33 d 3 . (22) 

Using the equations (21) and (22), and the fact that M = Midi + M 2 d 2 + M 3 d 3 , 
we obtain: 

Mi = J X 2 S 33 dS , (23) 
s 

M 2 = J-XiSndS, (24) 

s 

M 3 = j(X l S 32 -X l S 31 )dS . (25) 

s 

Finally, using the solutions for u-i given by the equations (19), and by the solution 
u 3 = (we are considering circular cross section so that (p(X 1 ,X 2 ) = 0) we 
can obtain the components of the strain tensor using the equations (8) and (17). 
Substituting the components of the strain tensor in the equations (7) and (16) we 
can obtain the expressions for the components S^i, S32 and S 33 of the stress tensor: 

S 31 = - f i(r)(k 3 -kf)X 2 , (26) 

S 32 =fi(r)(k 3 -k { 3 0) )X 1 , (27) 

S33 = E(r)((h - kf ] )X 2 - (k 2 - fc^Xx) . (28) 

The constitutive relations for the components of the total moment of the cross sec- 
tion can be written in a final form as: 

h 

Mi = (A* - fcfV J E(r)r 3 dr , (29) 
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h 




(30) 



o 



h 




(3D 



o 



where h is the cross section radius. 

Note that there is no constraint on the variation of the Young's or shear modulus 
with r (provided that the Poisson's ratio is constant). 

The constitutive relations given by equations (29-31) have the same form of the ho- 
mogeneous case. It means that within the approximations of the Kirchhoff 's theory, 
the static and dynamics of a rod are not affected by radial nonhomogeneities in its 
cross sections. Our calculations constitute a demonstration that a thin nonhomoge- 
neous rod behaves like a homogeneous one if the deformations have high radius of 
curvature as compared to the cross section radius. 

The equations (29-31) can be used to compare the rigidity of homogeneous and 
nonhomogeneous rods with the same cross section radius. By comparing the ex- 
pressions for Mi, M 2 and M 3 of homogeneous and nonhomogeneous cases, we 
can derive expressions for an effective Young's and shear moduli, E e j and fi e f, in 
terms of E(r) and /i(r), respectively. By effective we mean the values for Young's 
and shear moduli which result in the same values for M; (i = 1,2,3), if the rod was 
homogeneous. They are given by: 




(32) 



o 




(33) 



o 



11 



Consider the simple example of a rod like a coaxial cable, with the Young's modu- 
lus given by 



E(r) 



E for < r < r , 

(34) 

Ei for r < r < h , 



so that E and E 1 are the Young's moduli of the inner and outer parts of the cross 
section defined by the regions from the origin to r and from r to h, respectively. 
Substituting the equation (34) in the equation (32), we obtain the following expres- 
sion for the effective Young's modulus: 



i 

hj ' -M- h A 



Eef = E ( r 4) +25 1 [l-£) . (35) 



This simple expression can be used in various experimental set ups to find out 



the Young's m odulus of several systems. For example, consider the experiment of 



Salvadori et al. 



(2003a.b) where the Young's modulus of gold thin films deposited 
in Atomic Force Microscopy (AFM) cantilevers were measured. We cannot apply 
our method directly to their experiments because the cross sections of the can- 
tilevers are not circular, but our model provides another form of measuring the 
Young's modulus of gold thin films. If a thin film of gold (or another material) is 
made to grow on the cylindrical surface of a homogeneous rod, which Young's 
modulus E is known, by measuring the bending coefficient of the coated rod 
(which is related to E e f), the film thickness, h — r , and the total radius of the 
cross section, h, we can use equation (35) to obtain the Young's modulus, Ex, of 
the thin film. The only requirement is that the dimensions of the rod must be in 
accordance with the approximations of the Kirchhoff model. 

Another interesting application is the study of the rigidity of superficial layers of 
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nitride and carbonitride compounds within steel due to absortion of nitrogen by a 



process called nitriding (M iola et al. 
cations, in this case. 



1998). Our method can be used in two appli- 



In the first application, as for the case of rods coated by thin films, we can measure 
the Young's modulus of the substance of the layer of nitride and carbonitride com- 
pounds by producing rods of steel (which Young's modulus is known), submitting 
the rod to the nitriding process and measuring the bending coefficient of the ni- 
trided rod and the thickness of the layer. As in the previous case, the equation (35) 
gives the Young's modulus of the layer. 

The second application of our method to nitrided rods is an indirect measurement 
of the thickness of the layer of nitride and carbonitride compounds through the 
measurement of the Young's modulus of the nitrided rod. In this case, the Young's 
modulus of the steel, and of the substance that composes the layer, must be previ- 
ously known. We can obtain r through equation (35). 



In the Introduction, we mentioned that the Kirchhoff rod model has been used to 



study the elastic b ehavior of long pieces of DNA (Sc hlick . 



Tobias et al. 



1995 



Olson, 



1996 



2000). Since it is known that DNA is not a polymer with uniform 
cross section, our calculations guarantee that those studies about elasticity of long 
DNA's are not inconsistent with the use of the Kirchhoff rod model. Nevertheless, 
our approach could be also used to obtain information about the stiffness of the 
parts of the DNA. It is known that the most rigid part o f the double helix DNA is the 



phosphate backbone chain ( Calladine and D rew 



1999). The bases are connected, 



in a base-pair, through hydrogen bonds that are weak connections. Therefore, a 
complete base-pair in a DNA molecule is a net situation approximately close to the 
cross section of a coaxial cable. The inner part of DNA is formed by the bases that 
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are connected by weak hydrogen bonds and the outer part is formed by the strong 
phosphate backbone chains. If, by means of molecular approaches, it is possible to 
estimate the stiffness of the inner (outer) part of the DNA molecule, then the stifness 
of the outer (inner) part could be obtained, and the results could b e checked by 



expe rimentally measured stiffness of long DNAs (see, for example, (iSmith et al 

mm . 
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